The Current State of the Protected Apis mellifera mellifera Population in Russia: Hybridization and Nosematosis.

The Southern Urals of Russia are the habitat of one of the surviving populations of the dark forest bee-the Burzyan population of Apis mellifera mellifera. In this study, we present the results of the subspecies identification of bee colonies in the Altyn-Solok Nature Reserve in the Southern Ural Mountains using the intergenic mtDNA COI-COII locus and the assessment of the prevalence of nosematosis. Analysis of the mtDNA COI-COII intergenic locus in the studied sample showed that 30.4% of the colonies belong to the lineage C. The PCR diagnostics of nosematosis in 92 colonies selected from different sectors of the Altyn-Solok Nature Reserve showed that about half of the analyzed colonies were infected with Nosema apis. Nosema ceranae was found in eight colonies. Both of these factors can lead to the extinction of this population of the dark forest bee.

Improved Apis mellifera reference genome based on the alternative long-read-based assemblies.

Apis mellifera L., the western honey bee is a major crop pollinator that plays a key role in beekeeping and serves as an important model organism in social behavior studies. Recent efforts have improved on the quality of the honey bee reference genome and developed a chromosome-level assembly of 16 chromosomes, two of which are gapless. However, the rest suffer from 51 gaps, 160 unplaced/unlocalized scaffolds, and the lack of 2 distal telomeres. The gaps are located at the hard-to-assemble extended highly repetitive chromosomal regions that may contain functional genomic elements. Here, we use de novo re-assemblies from the most recent reference genome Amel_HAv_3.1 raw reads and other long-read-based assemblies (INRA_AMelMel_1.0, ASM1384120v1, and ASM1384124v1) of the honey bee genome to resolve 13 gaps, five unplaced/unlocalized scaffolds and, the lacking telomeres of the Amel_HAv_3.1. The total length of the resolved gaps is 848,747 bp. The accuracy of the corrected assembly was validated by mapping PacBio reads and performing gene annotation assessment. Comparative analysis suggests that the PacBio-reads-based assemblies of the honey bee genomes failed in the same highly repetitive extended regions of the chromosomes, especially on chromosome 10. To fully resolve these extended repetitive regions, further work using ultra-long Nanopore sequencing would be needed. Our updated assembly facilitates more accurate reference-guided scaffolding and marker/sequence mapping in honey bee genomics studies.

A revision of subspecies structure of western honey bee Apis mellifera.

The taxonomy of honey bee A. mellifera contains a lot of issues due to the specificity of population structure, features of biology and resolutions of honey bee subspecies discrimination methods. There are a lot of transition zones between ranges of subspecies which led to the gradual changes of characteristics among neighbor subspecies. The modern taxonomic pattern of honey bee Apis mellifera is given in this paper. Thirty-three distinct honey bee subspecies are distributed across all Africa (11 subspecies), Western Asia and the Middle East (9 subspecies), and Europe (13 subspecies). All honey bee subspecies are subdivided into 5 evolutionary lineages: lineage A (10 subspecies) and its sublineage Z (3 subspecies), lineage M (3 subspecies), lineage C (10 subspecies), lineage O (3 subspecies), lineage Y (1 subspecies), lineage C or O (3 subspecies).

Estimation of C-derived introgression into A. m. mellifera colonies in the Russian Urals using microsatellite genotyping.

BACKGROUND: Marker-assisted selection is well established in animal breeding method of selecting individuals with desirable traits in a breeding scheme based on DNA molecular marker patterns. OBJECTIVE: Genetic diversity and C-derived admixture into local purebred gene pool of A. m. mellifera colonies was assessed using polymorphism of nine microsatellite loci in order to provide further marker-assisted selection of desired honey bee colonies. METHODS: The genetic diversity and the level of C-derived introgression into A. m. mellifera colonies in the Shulgan-Tash Nature Reserve (Russia) was assessed based on nine microsatellite loci (ap243, 4a110, A24, A8, A43, A113, A88, Ap049, A28), which were analized using the fragment analysis of the PCR products in Applied Biosystems 3130 DNA Analyzer. Phylogenetic relationship of colonies was evaluated using Neighbor-Joining methods with Cavalli-Sforza and Edwards genetic distance using the PHYLIP 3.68. The model-based Bayesian clustering algorithm implemented in STRUCTURE 2.3.3 was employed to infer membership and introgression proportions (Q-value). RESULTS: In the Shulgan-Tash Nature Reserve colonies of A. m. mellifera subdivided into four groups by level of C-derived introgression. Only five colonies of A. m. mellifera had C-derived introgression which varied from 0.5 to 2%. The genetic diversity in colonies of A. m. mellifera varied from 0.12 to 0.40. The Neighbor-Joining tree demonstrates the genetic relationship of A. m. mellifera colonies, which subdivided into three groups with different levels of C-derived introgression. Group 1 combined five honey bee colonies Bort_1, Bort_2, Bort_3, Baisalyan_1, and Kush_7 with a fraction of introgression close to 0.000 and genetic diversity from 0.20 to 0.25. CONCLUSION: The results showed the excellence of nine microsatellite loci genotyping in estimation of genetic diversity, distinguishing the two European evolutionary lineages M and C and estimating C-derived introgression. These genetic parameters can be applied further to perform the marker-assisted selection of purebred dark European honey bees.